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Abstract 

Quantum statistical correlations and momentum distributions are calculated 
for a spherically symmetric, three-dimensionally expanding finite fireballs, 
for non-relativistic expansions applying plane-wave approximation. The new 
concepts of the geometrical temperature as well as the thermal radius are 
presented for the simplest case. The symmetry properties of the correlation 
function are shown to reflect the symmetry properties of the emission function. 
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Intensity interferometry has recently become a popular tool to infer the freeze-out surface 
heavy ion collisions and to study the space-time characteristics of the strong interactions 
in elementary particle reactions. The method was originally invented [|TJ to measure the 
angular diameter of distant stars or galaxies, the investigated objects being approximately 
static and the length-scales astronomical. In principle the same method is applied to measure 
the space-time characteristics of the strong interaction, where the objects under study are 
expanding systems, with life-times of a few fm/c (10 -23 sec) and length-scales of a few fm 
(10~ 15 m). 

In recent publications |||3| the new concept of the geometrical temperature and also 
the concept of thermally generated length-scales j|j5||||!| were introduced for relativistically 
expanding systems both in one and three dimensions. In ref. application for recent high 
energy heavy ion data from the CERN energy region [^,[7] was found: for three-dimensionally 
expanding, cylindrically symmetric, finite systems the Bose-Einstein Correlation Function 
(BECF) was shown to measure the thermal length scale, the region in the coordinate-space 
from where bosons with similar momenta may emerge. The size of this region is deter- 
mined by the freeze-out temperature, the freeze-out time and the flow gradients. It was 
shown that this thermal length-scale dominates the correlation function if the geometrical 
characteristic length scales are sufficiently large. However, the mt dependence of the pa- 
rameters of the BECF were shown to contain correction terms proportional to the ratio of 
the thermal/geometrical length scales. As a complementary effect, the momentum distribu- 
tion was shown to be dominated by the geometrical and not the thermal parameters, which 
provides us a possibility to measure the (large) geometrical scales by a systematic study of 
the single-particle momentum distribution. These effects appear for thermal lengths small 
compared to the geometrical size of the expanding system. The standard results |8||| are 
reproduced for the case the geometrical length-scales being smaller than the thermal ones: in 
this conventional case the BECF measures the geometrical size properly and the freeze-out 
temperature can be extracted from the momentum distribution. 

The motivation for the present Letter is to investigate how the thermal lengths and 
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the geometrical temperature reveal themself for the simplest systems where they appear at 
all. Since the thermal length-scales appear due the interplay of the flow-gradient and the 
temperature, a non-relativistically expanding thermal source with constant flow-gradient 
is the simplest possible case. The key feature is that the Bose-Einstein or Fermi-Dirac 
correlations appear due to an enhanced/decreased probability of finding identical particles 
in a similar momentum-state; i. e. the correlations vanish for large momentum difference. 
In case of an expanding source with flow gradients, the mean momentum of the emitted 
particles changes with changing locations. Locally, the width of the momentum distribution 
of the emitted particles is determined by the temperature; thus there will be a region in 
the coordinate space, where the change in the mean momentum of the emitted particles is 
of the same magnitude as the width of the local momentum distribution. The size of this 
region is referred to as the thermal length. If it is smaller than the geometrical size of the 
expanding system, the quantum statistical correlations (Bose-Einstein or Fermi-Dirac) shall 
be sensitive only to this thermal size, since this is the size of the coordinate-space from which 
particles with similar momenta may emerge. 

Thus the simplest case is a non-relativistic expansion with constant gradient. We may 
choose the case of the three-dimensionally expanding, finite systems with spherical symme- 
try, which is most probably realized in low energy heavy ion collisions. For a recent review 
on the intensity interferometry of low and intermediate energy heavy ion collisions (E/A < 
1 GeV) see ref. |0| . 



The particle emission is characterized by the emission function S(x,p) which is the 
probability that a particle is produced at a given x = (t, f) = (t, x, y, z) point in space-time 
with the four-momentum p = (E,p) = (E,p x ,p y ,p z ), where the particle is on mass shell, 
m? = E 2 — p 2 . The quantum- mechanical analogy to the classical emission function is the 



time-derivative of the Wigner-function, |11 



In terms of the time derivative of the Wigner-function both the momentum spectra and 
the quantum statistical correlation functions (QSCF-s) are prescribed. A useful auxiliary 
function is the Fourier-transformed emission function 



S(Ak, K) = J d A x S(x, K) exp(zA£; ■ x), (1) 

where 

Pi +P2 

Ak — pi — P2, and K = (2) 

2 

and Ak-x stands for the inner product of the four- vectors. Then the momentum distribution 
of the emitted particles, f{p) is given by 

f(p) = S(Ak = 0,K = p) J (3) 

which is normalized to unity, 

Vp/(p) = l. (4) 



The two-particle Bose-Einstein or Fermi-Dirac correlation functions (BECF-s or FDCF-s) 
are prescribed in terms of our auxiliary function, 

I S{Ak,K) | 2 

C(K,Ak) = 1 ± , (5) 

I 3{0, K) | 2 

where the + sign stands for bosons and the — sign for fermions, as was presented e.g. in 
ref. [fLl,9||. In this Letter the effect of final state Coulomb and Yukawa interactions shall 



be neglected as implicitly assumed by the above equation, and the time-derivative of the 
Wigner-function shall be approximated by classical emission functions. We assume com- 
pletely chaotic, thermal emission. For Fermi-Dirac correlations, the plane-wave approxima- 
tion is rather crude especially due to the attractive final state interactions in pp and nn 
systems for which systems a large data set was taken at low and intermediate energies ]TU 



The Coulomb repulsion in the pp channel plays also an important role. These final state 
interactions can be taken into account in the Wigner-function formalism along the lines of 
ref. |Tl|]. The Fourier-transformations which appear in this Letter are the consequence of 
the plane- wave approximation for the relative wave-function ipAki^yt). They correspond to 



integrals over | ipAk{r,t) | 2 = l + cos(AA; -x). In the case the final state interactions are taken 
into account, the Fourier-transformations are to be replaced by integral transformations over 
the 4>Ak{r, t) relative wave function which can be calculated for the Coulomb + Yukawa final 



state interactions as mentioned e.g. in [11 



For central heavy ion collisions at low or intermediate energies the target and the pro- 
jectile form a collective state (stopping) which can be described as a non-relativistically 
expanding fluid within the framework of hydrodynamical models. Due to the expansion, 
the fluid cools and at a certain temperature it freezes out. In case of full stopping and 
thermalization (spherically symmetric momentum distribution) the information about the 
initial direction is lost thus the final freeze-out density distribution becomes approximately 
spherically symmetric, too. 

We shall assume that the emission function is characterized by a given a distribution of 
production points J(r ) and by a distribution of the freeze-out times, H(t). The correlations 
between space-time and momentum-space shall be introduced by a non-relativistic momen- 
tum distribution. We assume that the expanding system at the freeze-out is rare enough 
so that the quantum-statistical single-particle distribution can be well approximated by a 
Boltzmann-distribution, 



f B {x;p) = Cexp 



^ (p — mu(x)) 2 ^ 
y 2mT j 



(6) 



Here u(x) is the (non-relativistic) velocity, the freeze-out temperature is denoted by T and 
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the quantity C = is a normalization constant. Thus the emission function is 

(2nmT) 3/2 

characterized as 

S(x;K) = f B (x;K)I(f)H(t). (7) 

In order to simplify the results we keep only the mean value and the width of the source 
distributions i.e. we shall apply Gaussian approximations for the distribution functions of t 
and r as follows 



f r 2 \ 



I(f 



H(t) 



exp 



(2nR 



1 



2R, 
t 



exp 



(t-t ) 



(8) 



(9) 



(2vrAt 2 ) 1 / 2 ^ 2 At 2 J 

We shall also briefly discuss our results for the general case. 

We select a velocity of the 3D expanding matter at space-time point x so that it be 
spherically symmetric and describe an expansion in all three directions with a constant 
gradient. Thus the velocity around the mean freeze-out time t is assumed to have the form 



u\x) 



(10) 



which describes a scaling solution of the non-relativistic hydrodynamical equations at the 



mean freeze-out time to f° r r << re f- [ 12fl - Further we assume that the freeze-out 
temperature in this non-relativistic case is smaller than the mass of the particles, T < m. 

With the help of the above flow pattern, the Boltzmann factor at the mean freeze-out 
time can be rewritten as 



J(r-r a )=f B (t ,r;K) = Cexp 



1 (r-r a {K)f\ 



\ 



2R% 



(11) 



i.e. the emission of the particles with momentum K is characterized by a Gaussian. The 
mean emission point f a (K) is proportional to the observed momentum and the width of the 
emission points around the mean is determined by the mean freeze-out time (inverse of the 
velocity gradient) and the freeze-out temperature: 



f a (K) = -K, 



m 



(12) 



This width Rt shall be referred to as the 'thermal radius'. 

When these approximations are valid, our auxiliary function can be rewritten as 
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r / - n K 

S(Ak, K) = H(AE) / d 3 f exp ( -i Ak f ) /(f) J(r - t„— ), (13) 

m 

where H(AE) stands for the Fourier-transformed freeze-out time distribution. Let's intro- 

— * 

duce the effective, K dependent particle emission function 

J^/(f) = /(f) J(f-t -). (14) 

m 

Thus for arbitrary, Fourier-transformable distributions the momentum distribution and the 
QSCF is given as 

d 3 n , 
oc f(p) = / d*rlf(r) = (15) 



dp x dp y dp z 

P 



J d?r /(f) J(f-t —), (16) 



m 



- ~ - I I e K ff (^k) | 2 
C(K, Ak) = 1± | H(AE) | 2 (17) 

|/^(0)| 2 

Thus the momentum distribution is given by a convolution and the QSCF shall be a sym- 
metric function in the CMS of the pair, where AE = 0. 

These general results take particularly simple form in the introduced Gaussian approxi- 
mations. The auxiliary function S(Ak, K) contains a factor which is the Fourier-transform 
of the product of Gaussians and so easily integrable. The two Gaussian factors can be 
combined into a single Gaussian by introducing the following notation: 



(r - r a ) r (r - r b ) + r c 

+ — = , (18) 

p2 p2 p2 

and the new variables can be determined by the requirement that the coefficients of the 
powers of f should be equal on both sides. This yields the following relations: 
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1 1 
— + — 



(19) 



IX. tXrp IX. 



G 



n 



P 2 P 2 



-> 2 ^2,^2 

r a n + r c 
R% 



(20) 



(21) 



Ri 



and now the factor exp 



( ?A 

1 r. 



V 



2R: 



becomes independent of the integration. This factor shall 



enter the momentum distribution and can be rewritten as 

( 

exp 



V 



2mT* 



(22) 



where the effective temperature T* is the sum of the freeze-out temperature and the "geo- 
metrical temperature" , 



T, = T + T G , where T G = m- 



R G 



(23) 



'-n 



The last line defines the geometrical temperature, being just the mass multiplied by the 

Rg 

square of the mean expansion velocity — . The geometrical temperature T G is related to 

to 

the geometrical radius R G exactly the same way as the freeze-out temperature T is related 
to the thermal radius Rt- 



d2 _ 2 Tg 

Rq = ^0 ' 



(24) 



m 



(25) 



m 



The remaining integral in the auxiliary function is trivial, resulting in the following 
momentum-distribution and correlation function: 



m 

C{K,Ak) 



(2irmT* 



,3/2 



exp 



\ 



2mT i , 



1 ± exp(-i? 2 AP) exp(-At 2 A£ 2 



(26) 
(27) 



The effective temperature T* and the effective radius parameter, i?* were defined in 



eqs. 



The relative momentum Ak, which enters the correlation function, is invariant under 
Galilei-transformations. However, the energy difference is not invariant even under the non- 
relativistic Galilei transformations. Let us re-formulate AE in such a way that its specific 
directional dependence be more transparent. The energy difference could be identically 
rewritten as 



AE 



Pi 



P2' 



' Pi+P2\ Ak 



2m 



= v K Ak, 

m 



(28) 



where we have introduced the mean velocity of the pair, vk 



P1+P2 K 



2m m 

Let us define the out direction to be parallel to the mean velocity of the pair, vk, and 
the perp index shall be used to index the remaining two principal directions, both being 
perpendicular to the out direction. This naming convention corresponds to the one used in 
high energy heavy ion collisions fll3",5fl. 

By this definition we evaluate the correlation function in a frame which is assigned to 
a given mean momentum of the particle pair and not that of the fireball: by changing the 
mean momentum we change our reference frame too. In this co-moving frame the relative 
momentum is decomposed as 



Ak i^Qouti QperpJ 



(29) 



and the correlation function can be rewritten as 



C (Q out jQ perp) 1 i exp ^ RperpQperp ^outQ 



2 

out I ' 



(30) 



Rperp — — > (31) 

R 2 out = Rl + v K 2 At 2 . (32) 

Thus the perpendicular components shall measure the effective radius R*. What is the 
interpretation of this quantity? 

There are two length-scales in the problem: the geometrical length-scale Rq and the 
thermal one Rt, the latter being generated by the flow-gradient and the temperature. The 



effective R* is dominated by the smaller of the two, cf. eq.|T9|. With other words, for 
large, relatively cold three-dimensionally expanding systems not the whole source can be 
seen by quantum statistical correlations, but only a part of the whole system, which is 
characterized by a thermal length-scale Rt- This corresponds to the limiting case Rt << Rg 
(i.e. T « Tq), resulting in a thermally dominated correlation function and a geometrically 
dominated momentum distribution: 

Rperp = Rti (33) 
R 2 out = R 2 T + v K 2 AT 2 , (34) 
T* = Tq. (35) 

The region, viewed by the QSCF is shown also on Fig.l. The generation of the effective, 
geometrical temperature for cold, expanding systems is illustrated on Fig. 2. 

It is interesting to investigate the other limiting case when R T » Rg- In this case we 
re-obtain the standard results: 

R lerp = R Gi (36) 

R 2 out = R 2 G + v K 2 Ar 2 , (37) 
T* = T, (38) 

i.e. if the thermal length scale is larger than the geometrical size in all three direction, the 
BECF measurement determines the geometrical sizes properly, and the momentum distri- 
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bution will be determined by the freeze-out temperature giving a thermal distribution for a 
static source. The region seen by the QSCF is shown on Fig. 3. 



From eq. ^ it follows that the out component in general shall be sensitive to the duration 
of the freeze-out time distribution, since it contains a term vx 2 At 2 . This term vanishes in 
the center of mass system (CMS) of the particle pair, since AE = vk = in CMS. The 
correlation function in the considered case becomes symmetric, when evaluated in the CMS 



of the pair, since 



Qi = Qler P + Qlut in CMS of the pair (39) 



is the invariant momentum difference. The correlation function shall take up the following, 
particularly simple form: 

C(Qj) = 1 ± exp(-RlQj) in CMS of the pair. (40) 

This is to be contrasted to the three dimensionally expanding, relativistic, cylindrically sym- 
metric case where the correlation function was found to be symmetric in the longitudinally 
comoving system (LCMS) of the particle pair ||, and not in their CMS. 

This result emphasizes the importance of the analysis of the symmetry properties of the 
correlation functions (CF-s), since in two cases we found that these properties reflect the 
symmetry properties of the emission pattern: an axially symmetric 3d emission generates a 
CF symmetric in LCMS while a spherically symmetric 3d emission function generates a CF 
symmetric in the CMS of the pair. In case the emission pattern describes a Id expansion 
the symmetry of the BECF is broken. Note that the symmetry appears in the 3d cases when 
the CF-s are dominated by the thermal length-scales and obviously a thermal distribution 
becomes symmetric if one can find the appropriately comoving reference frame. 

We have three input parameters: Rq, to and T which determine two measurables, R* and 
T*. Thus the geometrical parameters can be determined only as a function of one (within 
the model unknown) parameter, which we may choose to be T. 

In any case, bounds can be given since 
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> 4 Rl — 



ni 



(41) 



T(T* - T) 




Rl- > Rl 



(42) 



T 



Our model approximations are valid if m > T* which correspond to the requirement to > R*- 
In summary, we have calculated the momentum distribution and the quantum statistical 
(Bose-Einstein or Fermi-Dirac) correlation function for three-dimensionally expanding sys- 
tems in a non-relativistic model for spherically symmetric expansions, applying plane-wave 
approximation for a completely chaotic (thermal) source. 

In all the three principal directions, two length-scales are present: the geometrical and 
the thermal one. The latter is due to the combined effect of the flow-gradient and the 
temperature within the considered model for particle emission. The radius parameters of 
the correlation function are dominated by the shorter of the thermal and geometrical length 
scales. A complementary effect is that the effective temperature is a sum of the freeze-out 
temperature and the geometrical temperature. Thus the effective temperature is dominated 
by the higher of the two temperature scales. 

In the chosen case the Bose-Einstein (Fermi-Dirac) correlation functions were found to be 
symmetric in the center of mass system of the pair, had the geometrical sizes been larger than 
the thermal length-scale. Lower limits on the geometrical size and the freeze-out time can 
be obtained from a simultaneous analysis of the momentum distribution and the correlation 
function, eqs. [4^, [41|. 

We have given examples to show that an analysis of the symmetry properties of the quan- 
tum statistical correlation function might indicate the symmetry properties of the underlying 
flow pattern and density distribution. 
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Figure Captions 



Figure 1. Illustration of the region seen by the correlation function in case Rj* << Rg- 
The correlation function views the shaded area. In the general case the correlation 



function yields a radius i?*, given by eq. 19 



Figure 2. The generation of the effective temperature in case Rt « Rg- The particles 
at position r are emitted with a mean momentum K = mf/to. Since T << Tq in 
this case, the density distribution in space is mapped to a density distribution in 
momentum space. Thus the effective temperature shall be given by Tq = mR 2 G jt^. 

Figure 3. Illustration of the region seen by the correlation function in case Rt >> Rg- 
The correlation function views the shaded area. In the general case the correlation 



function yields a radius i?*, given by eq. 19 
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FIG. 1. 
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FIG. 2. 
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FIG. 3. 
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